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of J/ip and 4>(2S) transverse polarization at large transverse momentum is confronted with recent experimental 
data, and potential shortcomings of the theoretical analysis are discussed. 



1. Introduction 

The production of charmonium states at high- 
energy colliders has been the subject of consid- 
erable interest during the past few years. New 
results from pp, ep and e + e~ experiments have 
become available, some of which revealed dra- 
matic shortcomings of earlier quarkonium pro- 
duction models. In theory, progress has been 
made on the factorization between the short dis- 
tance physics of heavy-quark creation and the 
long-distance physics of bound state formation. 
The colour-singlet model has been superseded 
by a consistent and rigorous framework, based on 
the use of non-relativistic QCD (NRQCD) §], an 
effective held theory that includes the so-called 
colour-octet mechanisms. On the other hand, the 
colour evaporation model || of the early days of 
quarkonium physics has been revived B. How- 
ever, despite the recent theoretical and experi- 
mental developments the range of applicability 
of the different approaches is still subject to de- 
bate, as is the quantitative verification of factor- 
ization. Because the charmonium mass is still not 
very large with respect to the QCD scale, non- 
factorizable corrections || may not be suppressed 
enough and the expansions in NRQCD may not 
converge very well. 

In this situation cross checks between various 
processes, and predictions of observables such 
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as quarkonium polarization and differential cross 
sections, are crucial in order to assess the impor- 
tance of different quarkonium production mech- 
anisms, as well as the limitations of a particular 
theoretical approach. Among the specific predic- 
tions of NRQCD, transverse polarization of di- 
rect (i.e. not from B and \c decays) J/ip and 
tp(2S) hadroproduction at large transverse mo- 
mentum has emerged as the most distinct and 
most accessible signature ||. Transverse polar- 
ization also discriminates NRQCD from other ap- 
proaches, like the colour evaporation model which 
predicts the quarkonium to be produced unpolar- 
ized. In the following, I will briefly review the 
general mechanism of quarkonium production at 
hadron colliders in the NRQCD factorization ap- 
proach. The prediction of J/ip and ip{2S) trans- 
verse polarization at large transverse momentum 
is confronted with recent experimental data from 
the Tevatron, and potential shortcomings of the 
theoretical analysis are discussed. 

2. Quarkonium production at hadron col- 
liders 

In the NRQCD approach, the cross section for 
producing a quarkonium state H at a hadron col- 
lider is written as a sum of factorizable terms, 

Aa{pp -> H + X) 

= ]T da(pp -> QQ [n] + x) (O h [n]), (1) 

n 
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where n denotes the colour, spin and angular 
momentum state of an intermediate QQ pair. 
The short-distance cross section da can be cal- 
culated perturbatively in the strong coupling 
a s . The NRQCD matrix elements {0 H [n}) (see 
[|| for their definition) are related to the non- 
perturbative transition probabilities from the QQ 
state n into the quarkonium H. They scale with a 
definite power of the intrinsic heavy-quark veloc- 
ity v 0. (v 2 ~ 0.3 for charmonium and v 2 ~ 0.1 
for bottomonium.) The general expression ([!]) is 
thus a double expansion in powers of a s and v. 

The NRQCD formalism implies that so-called 
colour-octet processes associated with higher 
Fock state components of the quarkonium wave 
function must contribute to the cross section. 
Heavy quark pairs that are produced at short 
distances in a colour-octet state can evolve into 
a physical quarkonium through radiation of soft 
gluons at late times in the production process, 
when the quark pair has already expanded to the 
quarkonium size.0 

The production of S-wave charmonium in pp 
collisions at the Tevatron has attracted consider- 
able attention and has stimulated much of the 
recent theoretical development in quarkonium 
physics. The CDF collaboration has measured 
cross sections for the production of J /ip and 
t/>(2S) states not coming from B or radiative %c 
decays, for a wide range of transverse momenta 
5GeV < pt(ip) < 20GeV §,§. Surprisingly, the 
experimental cross sections were found to be or- 
ders of magnitudes larger than the theoretical 
expectation based on the leading-order colour- 
singlet model pcj . This result is particularly 
striking because the data extends out to large 
transverse momenta where the theoretical analy- 
sis is rather clean and non-factorizable corrections 
should be suppressed. The shortcoming of the 
colour-singlet model can be understood by exam- 
ining a typical Feynman diagram contributing to 
the leading-order parton cross section, Fig. 1(a). 
At large transverse momentum, the two internal 

2 Such a possibility is ignored in the colour-singlet model, 
where only those heavy quark pairs that are produced in 
the dominant Fock state (i.e. in a colour-singlet state and 
with the spin and angular momentum quantum numbers 
of the meson) are assumed to form a physical quarkonium. 



(a) leading-order colour-singlet: 
3 + S^ccfSiVs 
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(b) colour-singlet fragmentation: 
g + g^ [ccf Si 1 '] +gg]+g 




(c) colour-octet fragmentation: 
g + g^ccpsf]+g 




(d) colour-octet fusion: 




Figure 1. Generic diagrams for J/ip production 
in hadron-hadron collisions via colour-singlet and 
colour-octet channels. 

quark propagators are off-shell by ~ p\ so that 
the parton differential cross section scales like 
da/dp 2 ~ 1/Pt, as indicated in the figure. On 
the other hand, when pt 3> 2m c the quarko- 
nium mass can be considered small and the in- 
clusive charmonium cross section is expected to 
scale like any other single-particle inclusive cross 
section ~ \jp\- The dominant production mecha- 



3 



nism for charmonium at sufficiently large pt must 
thus be via fragmentation |ff[| , the production 
of a parton with large pt which subsequently de- 
cays into charmonium and other partons. A typ- 
ical fragmentation contribution to colour-singlet 
J/ij) production is shown in Fig. 1(b). While the 
fragmentation contributions are of higher order in 
a s compared to the fusion process Fig. 1(a), they 
are enhanced by a power p 4 /(2m c ) 4 at large pt 
and can thus overtake the fusion contribution at 
p t ^> 2m c . When colour-singlet fragmentation 
is included, the p t dependence of the theoretical 
prediction is in agreement with the Tevatron data 
but the normalization is still underestimated by 
about an order of magnitude (l2| , indicating that 
an additional fragmentation contribution is still 
missing. It is now generally believed that gluon 
fragmentation into colour-octet 3 Si charm quark 
pairs Jig], as shown in Fig. 1(c), is the dominant 
source of J/tp and tp{2S) at large p t at the Teva- 
tron. The probability of forming a J /tp{tp(2S)) 
particle from a pointlike cc pair in a colour-octet 
3 Si state is given by the NRQCD matrix element 
pj/^{2S)) fsf ] ]) which is suppressed by u 4 rel- 
ative to the non-perturbative factor of the lead- 
ing colour-singlet term. However, this suppres- 
sion is more than compensated by the gain in two 
powers of a s /Ti in the short-distance cross sec- 
tion for producing colour-octet 3 Si charm quark 
pairs as compared to colour-singlet fragmenta- 
tion. At 0{v A ) in the velocity expansion, two 
additional colour-octet channels have to be in- 
cluded, Fig. 1(d), which do not have a fragmen- 
tation interpretation at order a 3 but which be- 
come significant at moderate p t ~ 2m c |Q . The 
importance of the 1 sjf S> and 3 Pj contributions 
cannot be estimated from naive power counting in 
a s and v alone, but rather follows from the dom- 
inance of t-channel gluon exchange, forbidden in 
the leading-order colour-singlet cross section. 

The different contributions to the J/tp trans- 
verse momentum distribution are compared to 
the CDF data Q in Fig. 2. As mentioned above, 
the colour-singlet model at lowest order in a s fails 
dramatically when confronted with the experi- 
mental results. When colour-singlet fragmenta- 
tion is included, the prediction increases by more 



than an order of magnitude at large pt, but it 
still falls below the data by a factor of ~ 30. The 
CDF results on charmonium production can be 
explained by including the leading colour-octet 
contributions and adjusting the unknown non- 
perturbative parameters to fit the data. Numer- 
ically one finds the non-perturbative matrix el- 
ements to be of e>(lCT 2 GeV 3 ), see e.g. HQ, 
perfectly consistent with the v 4 suppression ex- 
pected from the velocity scaling rules. Similar 
conclusions can be drawn for tp(2S) production 
at the Tevatron. 
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Figure 2. Colour-singlet and colour-octet contribu- 
tions to direct J/tp production in pp — » J/tp + X at 
the Tevatron together with experimental data from 
CDF j^]. Parameters and NRQCD matrix elements 
as specified in (lq] , 

3. J/tp and tp(2S) polarization at the Teva- 
tron 

The analysis of the CDF data alone, al- 
though very encouraging, does not strictly prove 
the phenomenological relevance of colour-octet 
contributions because free parameters have to be 
introduced to fit the data. However, if factoriza- 
tion holds the non-perturbative matrix elements 
extracted from the J/tp and ■0(25') cross sections 
are universal and can be used to make predic- 
tions for various processes and observables. The 
single most crucial test of the NRQCD approach 
to charmonium production at hadron colliders is 
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the analysis of ip(2S) polarization at large trans- 
verse momentum. Recall that at large pt, ip(2S) 
production should be dominated by gluon frag- 
mentation into a colour-octet 3 Si charm quark 
pair, Fig. 1(c). When pt 3> 2m c the fragmenting 
gluon is effectively on-shell and transverse. The 
intermediate cc pair in the colour-octet 3 Si state 
inherits the gluon's transverse polarization and 
so does the quarkonium, because the emission of 
soft gluons during hadronization does not flip the 
heavy quark spin at leading order in the veloc- 
ity expansion. Consequently, at large transverse 
momentum one should observe transversely po- 
larized ip(2S) §. 

The polarization can be measured through the 
angular distribution in the decay ip — > l + l~ , given 
by dr/dcos 9 oc 1 + a cos 2 #, where 9 denotes the 
angle between the lepton three-momentum in the 
if) rest frame and the ip three-momentum in the 
lab frame. Pure transverse polarization implies 
a = 1. Corrections to this asymptotic limit due 
to higher order fragmentation contributions have 
been estimated to be small J!?]]. The dominant 
source of depolarization comes from the colour- 
octet fusion diagrams, Fig. 1(d), which are impor- 
tant at moderate pt. Still, at 0(v 4 ) in the velocity 
expansion, the polar angle asymmetry a can be 
unambiguously calculated within NRQCD 18 ll| 



in terms of the non-perturbative matrix elements 
that have been determined from the unpolarized 
cross section. In Fig. 3 we display the theoreti- 
cal prediction [Q (shaded band) for a in ip(2S) 
production at the Tevatron as function of the 
ip(2S) transverse momentum, taking into account 

3p(8) 



the 1 S U 8 ' > and fi P^' fusion channels and higher 



order corrections to the fragmentation contribu- 
tions. No transverse polarization is expected at 
Pt ~ 5 GeV, but the angular distribution is pre- 
dicted to change drastically as pt increases. A 
first measurement from CDF does not sup- 
port this prediction, but the experimental errors 
are too large to draw definite conclusions. 

The analysis of J/ip polarization is complicated 
by the fact that the experimental data sample |20tl 
includes J/ip that have not been produced di- 
rectly but come from decays of Xc and ip(2S) 
mesons. The contribution from the radiative de- 



cays of the higher excited states decreases, but 
does not eliminate, the transverse J ftp polariza- 
tion at large pt ]2l|j . Again, the NRQCD factor- 
ization prediction is not supported by the exper- 
imental data, see Fig. 4. 
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Figure 3. Polar angle asymmetry a for direct ip(2S) 
production in pp — > ip(2S)(—> p, + pT) + X at the 
Tevatron as a function of p± compared to data 
from CDF |o). See Q for details and the choice 
of parameters. Also shown is the prediction for 
two extreme choices of NRQCD matrix elements: 
(^[i^W^pW]) = o (upper dashed line) and and 
^[i s (8)]) = o, (Cf P (8) ]) = 5 • 10" 2 GeV 3 , 
(0*[ 3 S[ S) ]} = 5 • 1CT 4 GeV 3 (lower dotted line), as 
discussed in the text. 
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Figure 4. Polar angle asymmetry a for prompt J/tp 
production (i.e. not from B-decay) in pp — > J/ip{— » 
/i + /i~) + X at the Tevatron as a function of pt com- 
pared to data from CDF From 
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4. Discussion 

The absence of transversely polarized J/ip and 
ip(2S) at large transverse momentum, if con- 
firmed at the Tevatron Run II, represents a chal- 
lenge for the application of NRQCD factoriza- 
tion to charmonium production. In the follow- 
ing, I will discuss some possible uncertainties and 
potential shortcomings of the present theoretical 
analysis: 

• Higher-order QCD effects in the short-distance 
cross section should be included to improve the 
accuracy of the theoretical prediction, in partic- 
ular in the intermediate pt region. While no at- 
tempt has been made to estimate the impact of 
higher-order corrections on the spin-dependence 
of the cross section, preliminary studies indi- 
cate large NLO effects for the unpolarized pt- 
distributions |22Hl5|| . These corrections could 
strongly affect the determination of the NRQCD 
matrix elements from the Tevatron data. A full 
NLO analysis is however needed before quanti- 
tative conclusions can be drawn. The impact 
of non-vanishing transverse momentum of the 
incoming partons on the determination of the 
NRQCD matrix elements has been analysed using 
a Gaussian smearing of the pt distribution p2[ , 
Monte Carlo event generators [^3| and, most 
recently, the fc t -factorization formalism p^ , p5[ . 
The actual size of the effect, however, turns out 
to be very different for the different approaches 
studied in the literature. 

• The uncertainty in the determination of 
NRQCD matrix elements translates into an un- 
certainty in the predicted yield of transversely 
polarized J/tp and ip(2S). The extraction of 
NRQCD matrix elements from the unpolarized 
cross section at the Tevatron is affected by large 
theoretical uncertainties which have not yet been 
fully quantified, see e.g. [Jl5|,[l6). In Fig. 3 we 
also show the polarization prediction for the ex- 
treme choices of vanishing (O^^S^, 3 Pq 8 ^}) ma- 
trix elements (upper dashed line) and a very small 
{0^[ 3 S[ 8 ^]} matrix element (lower dotted line, see 
Figure caption). The latter choice is motivated 
by the recent analyses of J/ip and ip(2S) produc- 
tion at the Tevatron in the fc t -factorization ap- 



proach [£4 25 which predict a strongly enhanced 
colour-singlet contribution and a substantially de- 
creased value of (0^[ 3 s[ s) ]) ||. A decreased 
value of the (0^[ 3 S^]} matrix element, rela- 
tive to (O^sf \ 3 P (8) ]), would delay the on- 
set of transverse polarization and reduce the dis- 
crepancy between NRQCD predictions and data. 
However, such extreme values for the NRQCD 
matrix elements would indicate a flaw in our un- 
derstanding of the velocity scaling rules |7j] which 
imply that all three matrix elements are of ap- 
proximately equal size 0(v 4 ) « 0(1(T 2 GeV 3 ). 

• Heavy quark spin-symmetry is violated by 
higher-order terms in the NQRCD lagrangian and 
longitudinal polarization can arise if the binding 
of the charm quark pair into J/ip and ip(2S) pro- 
ceeds through two chromo- magnetic dipole tran- 
sitions. According to the velocity scaling rules 
these transitions are suppressed by v 4 and the 
corresponding corrections to transverse polariza- 
tion should not exceed « 5%. Predictions of 
heavy quark spin-symmetry for radiative charmo- 
nium decays agree quite well with experimental 
data H and there is no obvious reason why spin- 
symmetry should fail when applied to charmo- 
nium production. Still, the estimate of the spin- 
symmetry breaking corrections to transverse po- 
larization relies on the conventional velocity scal- 
ing rules 0| which have not been firmly estab- 
lished for the case of charmonium p6j. Alter- 
native velocity scaling rules J27p§| ] imply a hier- 
archy of NRQCD matrix elements different from 
the standard counting and might indicate larger 
spin-symmetry violating corrections. 

• Finally, charmonium production mechanisms 
beyond NRQCD could be responsible for J/ip and 
ip(2S) depolarization. For example, the leading- 
twist formalism of NRQCD factorization does not 
include possible rescattering interactions between 
the intermediate heavy quark pair and a comov- 
ing colour field. It has been shown that such 
rescattering corrections could yield unpolarized 
quarkonium |p9| . The analysis p9|| , however, re- 
lies on several simplifying assumptions, and fur- 
ther work is needed to establish the importance of 
comover interactions for charmonium production 
at large pt ■ 
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5. Conclusion 

The absence of a transverse polarization in J/ip 
and ip(2S) hadroproduction at large pt repre- 
sents a serious challenge for the application of the 
NRQCD factorization approach to charmonium 
production. Several potential problems with the 
current theoretical analysis have been discussed 
but no conclusive picture has emerged yet. If 
no transverse polarization is observed with higher 
statistics at larger values of pt in Run II of the 
Tevatron, one might have to conclude that the 
charm quark mass is not large enough for a non- 
relativistic approach to work in all circumstances. 

Note added 

After submission of this proceedings contribution, 
charmonium polarization at the Tevatron has been 
addressed in the fcj-factorization approach The 
authors argue that the 1 Sq 8 ^ production channel could 
be the dominant source of J/ip and ij)(2S) in the ex- 
perimentally accessible region of transverse momen- 
tum, which implies that quarkonium be produced 
mainly unpolarized. 
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